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Abstract 

Putting together kinematical and dynamical analysis, a complete study of 
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on the kinematical side, while on the dynamical side, Heavy Quark Effective Theory 
(HQET) is applied for an accurate determination of the hadronic matrix elements 
between the baryons A^ and A. Emphasis is put on the major role of the A;, 
polarization for constructing T-odd observables and the standard — lo mixing 
has the benefit effect of amplifying the process of direct CP violation between 
Ab and A^ decays. 
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1 Introduction 



Huge statistics of beauty hadrons are expected to be produced at the CERN-LHC proton- 
proton collider starting in 2007. Obviously this will lead to a thorough study of discrete 
symmetries, C, P and T in 6-quark physics, in the framework of the Standard Model (SM) 
as well as beyond the Standard Model. It is also well known that the violation oi CP 
symmetry via the Cabibbo-Kobayashi-Maskawa (CKM) mechanism is one of the corner- 
stone of the Standard Model of particle physics in the electroweak sector. Thanks to the 
foreseen CERN-LHC collider, non-leptonic and leptonic 6-baryon decays may allow us to 
get informations about the CKM matrix elements, analysis of the C — P—T operators may 
be performed and different non-perturbative aspects of QCD may also be investigated. 

Looking for Time Reversal (TR) violation effects in 6-baryon decays can provide us a 
new field of research. Firstly, this can be seen as a complementary test oi CP violation 
by assuming the correctness of the CPT theorem. Secondly, this can also be a path to 
follow in order to search for processes beyond the Standard Model. On the theoretical 
side, most of the time, the Time Reversal violation effect comes from an additional term 
related to the physics beyond SM and added by hand in the QCD Lagrangian. On the 
experimental side, various observables which are T-odd under time reversal operations 
can be measured, so that A^-decay seems to be one of the most promising channel to 
reveal TR violation signal. 

In a previous letter quoted as JJ, a general formulation based on the M. Jacob- 
G.C. Wick-J.D. Jackson (JWJ) helicity formalism has been set for studying the decay 
process Af, — > AV(1"), where A — pix^ with V l^l^ oi V ^ tt+tt^. Emphasis was 
put on the importance of the initial A^ polarization as well as the correlations among the 
angular distributions of the final decay products. On the dynamical side, the Hadronic 
Matrix Elements (HME) appearing in the decay amplitude were computed, at the tree 
level approximation, in the framework of the factorization hypothesis for two-body non- 
leptonic weak decay of heavy quark. A confirmation in the validity of the factorization 
hypothesis was the agreement found between the theoretical estimate of the A^ h^J/ip 
decay width and its experimental measurement given by the Particle Data Group (PDG) 
average [2]. 

In our present work, calculations are performed in a more exhaustive and detailed 
manner. On the dynamical side, both tree and penguin diagrams are involved in the 
evaluation of the HME. In our case, the non-leptonic Af,-decay proceeding through the 
W^-loop involves either the 6 — >■ u or 6 — > s transitions, respectively. The Cabibbo- 
Kobayashi-Maskawa matrix elements, Vub,Vus,Vcb,Vcs,Vtb,Vtd and Vts representing the 
charged current couphngs between quark transitions take place in the tree and QCD, 
electroweak, QCD-electroweak penguin diagrams calculated for the decay A;, AV{1^). 
Other diagrams such as the annihilation and box diagrams have been neglected in our 
approach. A key point in the calculation of non-leptonic baryon decays is the calculation 
of hadronic transition form factors between baryons that will be derived by making use 
of Heavy Quark Effective Theory (HQET). 

On the kinematical side, helicity methods are employed with an emphasis put on the 
polarization of the intermediate resonances A and V = J/ip, p'^jU and their particular role 
in testing the fundamental symmetries like Parity and Time Reversal operations as well. 
It is worth to say that the present work mainly based on the cascade-type analysis is very 
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useful in analysing polarization properties and Time Reversal effects. The main advantage 
of the decay cascade-type holds in a treatment where every decay in the decay chain is 
performed in its respective rest frame. Finally, full Monte-Carlo simulations including all 
the kinematics and dynamics are performed according to the computational model. 

The reminder of the paper is organized as follows. In Section 2, we present the 
kinematical properties of — > AV{1^) decays where the M. Jacob, G.C Wick and J.D. 
Jackson helicity formalism is extensively applied. In this section, we also emphasize the 
physical importance of polarization, mainly that of the intermediate resonances, A and V, 
respectively. In Section 3, we derive the decays of the intermediate resonances previously 
mentioned: all the calculations of A ^ P^r^, V l^l^ and V ti^tt~ are performed a 
la Jacob- Wick-Jackson where the role of the helicity frame is strongly alighted. 

After the kinematical analysis of A^ — > KV{1') decays, we focus on the dynamical 
analysis of A^-decay in Section 4. Transition form factors are derived in Heavy Quark 
Effective Theory which is well suited for A^-decay since the h quark mass is large enough 
to play the adequate energy scale in order to make an expansion in l/mt of the QCD 
Lagrangian written in the HQET formalism. Corrections to the order of 0{l/mb) are 
included into the weak transition form factors between the baryons A^ and A and an 
usual model for describing baryon wave functions is also applied in our calculations. The 
weak decay amplitude of the analysed process is derived by following an Effective Field 
Theory (EFT) approach. Operator product expansion and Wilson coefficients are used 
in order to include long and short distance physics, respectively. Moreover, in the special 
case oi V — p*^ and cu, the charge symmetry violation mixing between these two vectors 
is included since it may rise to a large signal of CP violation due to a large strong phase 
dynamically produced at the ou resonance. 

In Section 5, we list all the numerical inputs, CKM matrix elements, quark and meson 
masses and decay constants that are needed for calculating all the physical observables 
related to our analysis. Section 6 is devoted to results and discussions regarding all 
the simulations we made for A;, — > AV{1~) decays: branching ratio, BTZ, CP violating 
asymmetry parameter, acp, helicity asymmetry parameter, aAS, polarization, P, and 
various angular distributions. In this section, we also stress on the main physical 
observables which can be measured with the future detectors at the CERN-LHC proton- 
proton collider. In Section 7, we discuss in detail the search for Time-odd observables and 
we conclude this section by proposing an observable which is odd under Time Reversal 
operation: the distributions in cosine and sine of the normal-vectors to the decay planes 
of the intermediate resonances, A and V, in the A5 frame. Finally, our main results are 
summarized in the last section where conclusions and global search for discrete symmetries 
are discussed. 

2 Kinematical properties of A5 — ^ AV{1~) decays 

The hyperons produced in proton-proton collisions as well as in other hadron collisions 
are usually polarized in the transverse direction. The average value of the hyperon spin 
being non equal to zero and, owing to Parity conservation in strong interaction, the spin 
direction is orthogonal to the production plane defined by the incident beam momentum, 
Pi, and the hyperon momentum, p^. Usually, the degree of polarization depends on the 
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center of mass energy, a/s, and the hyperon transverse momentum ISj. In the case of 
pp —>■ Af, + X ai high ^/s, the 6— quark could replace the s— quark of an ordinary A and 
it is expected that the beauty baryon, Af,, will be polarized in the same way than the A 
hyperon. We define Np as the normal vector to the production plane: 

iVp = A^. (1) 

bi X Pb\ 

where pi and pb are the vector-momenta of one incident proton beam and A;,, respectively. 

Let (AbXYZ) be the rest frame (see Fig. of the A^ particle. The quantization axis 
(AfoZ) is chosen to be parallel to Np. The other orthogonal axis (Af,X) and (A^y) are 
arbitrary in the production plane. In our further analysis, the (A(,X) axis is taken parallel 
to the momentum pi, which constrains the direction of the {AbY) axis in the production 
plane. The spin projection, M,, of the A;, along the transverse axis {AbZ) takes the values 
±1/2. An important physical parameter is the A^ Polarization Density Matrix (PDM), 
(p^'). It is a (2 X 2) hermitian matrix with real diagonal elements verifying X]i=i Pu' ~ ^■ 
The physical meaning of the diagonal matrix elements is the following: the probability of 
getting Ah produced with Mi = ±1/2 is given by p^i and P22, respectively. Finally, the 
initial Af, polarization, Va^, is given by (^a^ ■ Np) = Va^ = Pn — P22 • 

In the framework of the JWJ helicity formalism [4^, the decay amplitude, AQ{Mi), 
for Ab{Mi) —> A{Xi)V{X2) is obtained by applying the Wigner-Eckart theorem to the 
cS-matrix element: 

Ao(M,) = (l/2,M,|5(°)b,e,0; Ai, A2) = A^x.M^b ^ AV)D]/,%^{^,9,0) , (2) 

where p = {p,6,(f)) is the vector-momentum of the hyperon A in the A;, frame (Fig. P). 
Ai and A2 are the respective helicities of A and V with the possible values Ai = ±1/2 
and A2 = —1, 0, ±1. The momentum projection along the (A) axis (parallel to p) is given 
by Mf = Ai — A2 = ±1/2. The Mf values constrain those of Ai and A2 since, among 
six possible combinations, only four are physical. If Mf = ±1/2 then (Ai,A2) = (1/2,0) 
or (-1/2,-1). If Mf = -1/2 then (Ai,A2) = (1/2,1) or (-1/2,0). On other side, the 
hadronic matrix element, A(Xj^^x^){Ab — > AV), contains all the decay dynamics describing, 
for a set of helicities Ai and A2, the hadronic part of the decay A;, AV and finally the 
Wigner matrix element, 

fjs a n\ — 



'MiMf 



,e,0) = di,^,,^.ie)exp{-tMi(j)) , (3) 



is expressed according to the Jackson convention 4J. According to standard rules, the 
differential cross-section, da, must take account of the undetermined helicity initial state 
and include a summation over individual final helicity states, the total angular momentum 
along the helicity (A) axis, A = Mf = Ai — A2, being fixed. We get the expression: 



d^^J2 J2PM.Mi\Ax.M)i^b - AV)fd]!,yj,lexpz{M: - M,)0 . (4) 

As it is known. Parity is not conserved in weak interactions therefore the hadronic matrix 
element A(Xi,X2){M AV) is not equal to A(-\j^-\^){Ab AV). In order to handle in an 
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easy manner lengthy calculations, the following mathematical expressions are introduced: 
7(±l/2) = Yl PM,M^4ii/24^i/2expi(M; - M,)0 . (5) 

Mi,Ml 

Noticing that V^^ — p++ — p^^L and p^\. + pt^_ — 1 (i.e. normalization condition), the 
expression written in Eq. (5) become simplified: 

7(±1 /2) = ^ (l ± P^'-cos e ± 23fJe(p+^_ exp i(f))sm 9^ . (6) 
Other expressions are also introduced, like: 

(I;(+l/2) =|Ai/2,i)(A;, ^ AV^)|%(-l/2) + |Ai/2,o)(A6 ^ AV^)|%(+l/2) , 
u{-ll2) =|^(_i/2,-i)(A, ^ Ay)|%(+l/2) + |A-i/2,o)(A6 ^ Ay)|%(-l/2) , (7) 

which represent respectively the weight of the A helicity states +1/2 and —1/2 along the 
(A) axis. With such notations, the differential cross-section can be rewritten in the simple 
form: 

^ = AT (^(+1/2) +^(-1/2)) , (8) 

where A" is a normalization constant. It is worth to underline that, apart from the 
four hadronic matrix elements, A(\^^\^){h.i, — > AT^), which have a dynamical origin, 
the above cross-section needs three real parameters in order to be fully determined: 
7^^*, 3?e(p^'L) and '^m{p^_). The relation shown in Eq. (8) can be put in a more compact 
form by introducing the helicity asymmetry parameter, a^s, which is related to the final 
hehcity value Mf = A = Ai — A2. By means of the following relations, 

|A6(+)|' =|^(i/2,o)(A6 A\/)|2 + |^(_i/2,-i)(Ab ^ A\/)|2 , 

|A,(-)|2 =IA-i/2,o)(A;, ^ + IAiAi)(A6 ^ , (9) 

and defining the helicity asymmetry parameter, a as, ^ 

°^'=iM+,i^+iM-)r 

the final expression of the differential cross-section reads as 

^ ocl + aAsV^'cos e + 2aAs'^e{pX'_ exp i0)sin 9 . (11) 

Then, by averaging over the azimuthal angle, 0, a standard relation is obtained for the 
polar angular distribution: 

(xl + aAsV^''cos9 , (12) 



dcos9 



where it can be noticed that polar angular dissymetries are intimately related to the initial 
polarization of the Aj, resonance. 
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2.1 Physical importance of the polarization 



Parity violation in weak decays into necessarily leads to a polarization process of 
the two intermediate resonances A and V. In order to determine the vector-polarization 
of each resonance, a new set of axis is needed according to which properties of the vector- 
polarization V^^^ with i = A,V by Parity and Time-Reversal (TR) transformations will 
be more obvious. Let ez be the unit vector which is parallel to the preceding vector Np, 
Np being transverse to A;, production plane, and p the momentum of a resonance in the 
Afo rest frame. The following unit vectors are defined (the index (i) being dropped): 

P ^ ez X cl ^ ^ , . 

ei = - , ct = ^ , cn = cl X ct . (13) 

P \ez X clI 

In this new frame, the vector-polarization of any resonance defined in the original A^ 
frame can be rewritten like: 

^ p^W^^ + p^i^)g^ + p^Wg-^ ^ (14) 

such that each component can be easily computed, Pj^^ = P^*^ ■ ej with j = L, A^, T. 
We recall that Pl, Pn, Pt are named longitudinal, normal and transverse polarizations, 
respectively. Using the properties of the spin, s, and the vector-polarization, V, under 
Parity and TR operations, we can deduce those of Pj*^ displayed in Table [H 

It is worth noticing that the normal polarization (of any resonance) P^r is odd by Time- 
Reversal operation. It allows us to say that the knowledge of the vector-polarization, V, 
and particularly its normal component, P/v, are essential observables for performing any 
test of validity of Time-Reversal symmetry. 



2.2 Polarization of the intermediate resonances 

Intermediate resonance states, A and V, can be described by a density-matrix named 
which analytic expression is given by standard quantum-mechanical relations: 

pf = r^p^^T , with Tr{pf) = ^ = NW{e, 0) , (15) 

where T is the transition-matrix related to the iS-matrix by 5 = 1 + iT. Its elements 
being the hadronic elements, A(Xj^^x^){Ah AV), mentioned beforehand. With the help of 
the density- matrix p^ , the vector-polarization of the global system made out of A and V 
can be estimated according to the relation: 

P=(S) = ^'^^'^^ (16) 

where Tr is the trace operator and S the spin of the global system defined by = 15*1 + 521, 
5*1 and 5*2 being the spin of the baryon A and the vector meson V, respectively. 

The striking point is that the final system is composed of two correlated subsystems 
because only four spin-states instead of six are required to describe the quantum system 
{A (B V). So, in order to determine the kinematics features of each resonance, namely 
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A and V, and particularly there vector-polarization, and , some modifications 
of the previous relations are requested. Our task will consist in estimating individual 
density-matrix and vector-polarization for each resonance by strictly applying quantum- 
mechanical rules 13 El • 

2.3 Determination of 7^^ 

A helicity states and its vector-polarization can be extracted from the above relations by 
summing over the kinematics variables of the vector-meson states, especially its helicity 
states. Departing from formulas (15) and (16), the following relation can be inferred: 

V'^Wie, 0) = NJlxie, 0; X\p^a\e, 0; A) . (17) 

This clearly shows that A polarization depends a priori on its emission angles. In Eq. (fT7j) . 
we have also defined: 

• A = Ai — A2, the vector-meson helicity A2 which varies with the kinematical constraint 
X = Mf = ±1/2. 

• a are the Pauli spin 1/2 matrices and they are defined in the (Axyz) rest-frame where 
the axis x, y and z are respectively identified with the normal, transverse and longitudinal 
axis previously defined. 

The matrix elements given in the right-handed side of Eq. (|17|) can be explicitly 
calculated [7] and the three components of get the following expressions: 

PtW{e,(j)) oc cj(+l/2) - c<;(-l/2) , 

1^(^,0) cx 2g?e((^,0,l/2|p^|0,0,-l/2)) , 
P^W{e,(j)) oc -2S5m((^,0,l/2|p^|^,0,-l/2)) . (18) 

Explicit calculations lead to: 

PtW{e,<f)) oc 7(+l/2)(lAi/2,o)(Af.^A^)r-|A-i/2,-i)(A6^AF)f 

+ 7(-l/2)(|AiAi)(A.-AV^)P-|A-i/2,o)(A5-AF)r) , (19) 



PjO W{d, 0) oc (^AiA0)(Afe AV)Al_y,^,^{At ^ AV) 

(-V^'^ sine + 23fJe(exp (i0)p+'_) cos^ + 2i53m(exp (i0)p+'_)) ] , (20) 

and 

P^W^(e,0) (X -3m^^(i/2,o)(Afe^AV^)^^_i/2,o(Ab^AV^) 

(^-V^b sine + 23?e(exp (i0)p+'_) cos^ + 2i5>m(exp (20)p+'_)) j . (21) 
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In these detailed formulas given in Eqs. (19-21), we underline the importance of two input 
parameters: 

(i) the initial A;, polarization, V^'', and (ii) the non-diagonal matrix element, p^^L, which 
appears essentially in the interference terms. 

2.4 Longitudinal polarization of V 

Vector meson has spin S(^2) — 1 and therefore three hehcity states. The components of 

its vector-spin operator, Sj,j = x,y,z, are (3 x 3) matrices and they are defined in a 
new set of longitudinal, normal and transverse axis. The longitudinal one being the spin 
quantization axis and verifying the standard relation: 5*^1 A2) = A2IA2). 

Computation of the normal and transverse components, respectively and , as 
well as the analytical method of extracting these values from experimental data will not 
be exposed here. In this paper, in a first aim, only the longitudinal component, P^, will 
be estimated. According to the basic relation: 

V"" W{e, 0) = AT E;,, (Ea, {e, 0, Ai, A2|p^^|^, <P, Ai, A2)) , (22) 
can be deduced in a relatively easy way and one gets, 

W{0, <t>) = 7(-l/2)IAiAi)(A;. - m? - 7(1/2)1 A-iA-i)(A6 ^ (23) 

where the kinematics functions 7(±l/2) have been defined in Section 2. It may be noticed 
that the mathematical form of P^ is very similar to that of P^ (same longitudinal axis) 
so that the helicity value A2 = does not play any role in the longitudinal polarization 
of the vector-meson. 

3 Decays of the intermediate resonances 

Performing appropriate rotations and Lorentz boosts, the decay products of each 
resonance can be studied in its own helicity frame (see Fig. 1). The quantization axis 
of each frame is chosen parallel to the corresponding resonance momentum in the A^ 
frame i.e. Oi^iUpa and O2-22IIPV — —Pk- Again, conservation of the total angular 
momentum (in each frame) constrains the hclicitics of the final particles. For the 
decays A(Ai) P(A3)7r-(A4) and V{\2) ^ or V{\2) ^ {\^)h+ {X^^) , 

the respective helicities are (A3, A4) = (±1/2,0) and (A5, Ae) = (±1/2, ±1/2) in case of 
leptons or (A5, \q) — (0, 0) in case of 0~ mesons. 

In the A helicity frame, the projection of the total angular momentum, m^, along the 
proton momentum, pp, is given by mi = A3 — A4 = ±1/2. In the vector meson helicity 
frame, this projection is equal to m2 = A5 — Ae = — 1, 0, +1 if leptons and m2 = if pions. 
The decay amplitude, /Ij(Ai), of each resonance can be written similarly as in Eq. (2), 
requiring only that the kinematics of its decay products are fixed. Thus, we obtain, 

Ai(Ai) =(Ai, mi|^(^)|pi, ^1, 0i; A3, A4) = Aa3,A4)(A ^ p7r-)Dif„\(0i, ^1, 0) , 

^2(A2) =(A2,m2|5(')b2,^2>2; As, Ae) = A^x,m){V ^ /i+/i-)Di:^,(02, ^2, 0) , (24) 
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where 6i and 0i are respectively the polar and azimuthal angles of the proton momentum 
in the A rest frame while 62 and 02 are those of i~{h~) in the V rest frame. Finally, 
^(A3,A4)(A PTT") and A(Xg^xQ){V — > ,h'^h~) are the dynamical amplitudes of the 
strong interaction part of the decay processes A — > pn" and V — > h^h' , respectively. 

3.1 Analytical form of the decay probability 

The general decay amphtude, Aj{Mi, Ai, A2), for the process of sequence decays Aft(Mj) — > 
A{Xi)V {X2) p7r~£^i~ (h^hr) can be factorized according to the three decay amplitudes 
y4o(Mj), Ai(Ai) and A2{\2)- It must include all the possible intermediate states, so that 
a sum over the helicity states (Ai, A2) is performed: 

A%Mi, Ai, A2) = Yl AoiMi)Ai{X,)A2{X2) ■ (25) 

Ai,A2 

The decay probability, da, depending on the amplitude, A%Mi, Ai, A2), takes the form, 
da^Yl PM,Mi^KMi, Ai, A2)^f (M,, Ai, A2) , (26) 

Mi,Ml 

where the polarization density matrix, P^.mi, is used to take into account the unknown 
initial A;, spin states, Mi. Since the hclicitics of the final particles are not measured, 
a summation over the helicity values A3, A4, A5 and Ag is performed as well. In a final 
step, the decay probability, da, can be written in a such way that only the intermediate 
resonance helicities appear: 

da<x Yl D^.-x,,K-W,{0, 0)Aa„a.)(A(. ^ AF)^^,, ^,,)(A, ^ kV) 

^Aai(^l>l)GL^(^2,02), (27) 

where F^^^, {9i, 0i) and G\^^,^ [62, ^2) are the hadronic tensor elements describing the decay 
dynamics of the intermediate resonances A — > p7r~,V t^l''{h'^h~) and -Da,a' is a 
kinematical factor. Their expressions are: 

FtwS^iAi) = (lAiA0)(A-P7r-)|Xfv2(^i)4fi/2(^i)+ 

I A-i/2,o)(A ^ P7r-)|'dif_,/2(^i)<-i/2(^i)) expi(A; - AO0i , (28) 

and, 

Gl4^2.ct>2) = \Ax.M){y - ^^^",^^^-)l'<™.(^2)4>.(^2)expi(A', - A2)02 , 

-^5,-^6 

^A,-A.Ai-A^(^, ct>)^Y PM,M^4'A.-A.(^)<,Ai-Ai(^)exp^(M; - M,)0 , (29) 

with 777.2 = A5 — Ae. As far as the vector- meson V is concerned, in case of lepton pair 
in the final state and because of parity conservation, two hadronic matrix elements. 
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^{±i,±i)(y ^ l^l^ , h^h^), are necessary whereas only one, ^(o,o)(^ l^l^ , h^h^), 
is required in case of pseudo-scalar mesons. 

From the above relations, angular distributions of the decay products of both the two 
resonances can be deduced. Although these angular distributions are well established in 
the literature, the method outlined before allows us to compute the Polarization Density 
Matrix (PDM) elements of each resonance. We underline that these angular distributions 
are important ingredients for the determination of its vector-polarization, P*^*^ (^ = A,V). 



3.2 A — > piT decay 

Departing from Eq. (jTfjl . integrating over the angles 6,(j),62 and 02, and summing over 
the vector helicity states, the general formula for proton angular distributions, Wi{9i, 0i), 
in the A rest-frame can be obtained: 



(P++ + P--) + (P++ - Pt.)a% cose, - -V^-a\s^e 



sin6'i 



(30) 



where the PDM elements, p^ , of the A hyperon are not normalized yet and expressed as: 



I Gj;(^2,02) |^(±1/2,0)(A6^AV-)|'+ (31) 

\ Gli±i(^2,02)IA±iA±i)(Afe ^ , 

J d2,4'2 

G^oi02, 02)A-iAo)(A, -> Ar)^^i/2,o)(A6 ^ AV) . (32) 



d2,4>2 



Then, integrating over the azimuthal angle, 0i, the standard expression for the polar 
angular distribution can be recovered: 



da 



oc 1 + a^gV^cosOi 



d cos di 



(33) 



where is the usual A helicity asymmetry parameter and is the A polarization 
according to the helicity axis. Integrating the formula given in Eq. ()27p over the polar 
angle 6'i, it allows us to deduce the proton azimuthal ang'u/ar distribution: 



da 



TC 



OC 1 + —V 



P12 exp {i4>i 



(34) 



It will be shown later that the PDM element, p^g' is real in the framework of the 
factorization hypothesis used to compute the hadronic matrix element ^Ai,A2(A;) — > AV). 



3.3 V £+£-, h+h- decays 

Vector meson, V, decaying into a lepton pair or a hadronic one is described by the (3 x 3) 
hermitian matrix G^^ {62, 02)- The angular distributions, W^2(^25 ^2); in the V rest-frame. 
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are obtained by integrating Eq. ()27|) over the angles d,(j),6i and (pi. Summing over the 
two A hehcity states, one gets, 



1^2(^2,02) OC (p]^ + p;;.)(GS;(^2,02) + Gll±l(^2,02)) 



4 



-.V 



Pa exp [1^2) 



sin 2^2, (35) 



where the PDM elements, pjj, of the meson V are (to a normalization factor): 



V 

PiJ 



L 



5a2A^IA±1/2,0)(A6 ^ A^)|' + ^A2±A(,IA±l/2,±l)(Afe ^ A\/)|' 

1^(1/2,0) (Afe ^ AV)^^i/2,i)(A6 ^ AV) + /i.e.} 



A-1/2,0)(A6 ^ AV)^^_i/2,_i)(A;, AV) + /i.e. 



V^hh{ll) 



(36) 



where, J^v->hh{ii), takes the following form according to the given decay: 
Mv^hhiii) =|^(o,o)(^ ^ h+h-)]"^ , for 1/ ^ , 

Mv^hhiii) =IAiA-i/2)(^ ^ ^+^~)l' - 2|^(+i/2,+i/2)(F ^ /+/-)|2 , for K ^ /+/- . 

Taking advantage of parity conservation in strong or electromagnetic K-decays, a relatively 
simple polar angular distribution can be obtained after averaging over the azimuthal angle 
02- Thus, in case of leptonic decays: 

da 



OC 



d cos 02 

while, for hadronic decays, one gets: 

da 



(l-3po'o)cos2^2 + (l + Po'o) 



(SPm - l)cos^^2 + (1 - Poo) , 



(37) 



(3^ 



d cos 6*2 

where is the PDM element indicating the probability for the vector-meson V to be 
longitudinally polarized. In both the two angular distributions Wi and W2 displayed in 
Eqs. (30) and (35), it is worth noticing the major role played by the Ab polarization, V^'', 
especially in the azimuth angular 0i and 02 distributions, which arise from interference 
terms. 



4 Hadronization 

In this section, the Heavy Quark Effective Theory jHUHlEj (HQET) formalism will be used 
to evaluate the hadronic form factors involved in A^-decay. Weak transitions including 
heavy quarks can be safely described when the mass of a heavy quark is large enough 
compared to the QCD scale, Aqcd- Properties such as flavour and spin symmetries 
can be exploited in such way that corrections of the order of l/mg are systematically 
calculated within an effective field theory. Then, the hadronic amplitude of the weak 
decay is investigated by means of the effective Hamiltonian, AB = 1, where the Operator 
Product Expansion formalism separates the soft and hard regimes. 
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4.1 Transition form factors 



The decay, A;, AV, is usually described by the following general amplitude, A{Ab 
AV), 



A{Ab ^ AV) = iUa{pa, SA)el A^^-f^ + A^ p^-f^ + Vi-f^" + V2 f/A,(pA„ sk. 



(39) 



where, Vi and Aj, are respectively the vector and axial components of the transition form 
factors characterizing the decay A;, — > AV . The momentum and spin of the baryons A 
and Afo are given by pa^^j and sA(j,), respectively, e*^ is the vector polarization and U^f, and 
Oa are the Dirac spinors. Based on Lorentz decomposition, the hadronic matrix element, 
(A(]9a, sa)|s7^(1 - 75)6|Afe(pA6, saJ), may also be written as, 

(A(pa,sa)|s7^(1 -75)&|Af,(pA6,SAj) = 

t/A(PA, sa) [(/i(g')7/. + ^/2(g')cT;..g" + /3(g')g^) 

- (^i(9^)7m + W2{q^)(Tf,uq'' + 93{q'^)qp.)l5 Ua^{pa^, saJ , (40) 

where, = {pa^ — Pivfi defines the momentum transfer in the hadronic transition. 
The form factors, fi{q^) and gi{q'^), refer to the vector and axial parts of the transition, 
respectively. 

Another way of parameterizing the electroweak amplitude in decays of baryons is the 
following: 



(A(pA, sa)\sY{^ - l5)b\Ab{pA„ saJ) = 

f/A(PA, Sa) \(F,{q')r + F2{q')vl + Fs{q')^) 



Gi{q'h^ + G2{q')vl + Gs{q')^) 



75 



UaSpa^sa,). (41) 



By comparing the two sets of form factors given in Eqs. (pn|) and ()41|) . one gets the 
relations between the /i(g^)'s {gi^q'^Ys) and Fj(g^)'s (Gj(g^)'s) as follows. 



Mq') 
f2{q') 

Mq') 



Fi(g^) + (mAi, + rriA) 
+ 



2mA. 2mA 



F2{q') , F,{q^ 



2mA, 
F2{q^) 
2mA, 



2mA 

F?Xq') 

2mA 



(42) 



and. 



9i{q'^) = Gi{q^) - (mA, - mA) 



G2{q') , G,{q^) 
2mA, 2mA 



92[ 



G2if)_^Chif) 
2mA, 
G2{q') 



2mA, 



2mA 

g3(g^) 

2mA 



(43) 
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In case of working in the HQET formalism, the matrix element of the weak transition, 
Af, A, takes the following form, 

(A(PA, sa)|s7^(1 - 75)6|Afe(pA,, saJ) = Ua{pa, sa) k(g') + KHq'')] Ua,{pa„ saJ • (44) 



In Eq. (|44|) . va^, defines the velocity of the baryon A^. Writing the momentum, p^^^, of 
the heavy baryon. A;,, such as, 

PAt = rubVA^ + k , (45) 

where, k, is the residual momentum, the velocity, va^, of the heavy quark, b, is almost that 
of the heavy baryon, A5. Because of rrih ^ Aqcd, the parameterization of the hadronic 
amplitude, A{A{, AV), in term of the velocity, v^^, gives us a reasonable picture where 
corrections only arise in the l/rrih expansion. 

The Stech's approach ^T] mainly assumes that the off-shell energy of a constituent 
diquark is close to its constituent mass (without dependence of its space momentum) in 
the rest frame of the baryon. Moreover, the spectator quark retains its original momentum 
and spin before final hadronization. Therefore, the energy carried by the spectator quark 
is equal to that of the spectator in the rest frame of the final state particle and the 
relevant 6-quark space momenta are much smaller than the b quark mass: indeed, it 
is assumed to be of the order of the confinement scale, Aqcd- This approach firstly 
used in the meson case can be generalized to a heavy baryon considered as a bound 
state of a 6 quark and a scalar diquark. Thus, in the baryon case, the hadronic matrix, 
(A(pa, '5a)|s7^(1 — 75)6|Ab(pAb5 saJ), written in terms of components of Dirac spinors such 
as, 

Us{PA,ms)'^^,{l - 'j5)Ub{pAb = 0,mfe) , 

leads to the following expressions for the form factors, 9i{q^) and ^2(0'^), when the heavy 
quark mass goes to infinity: 

9i{q^) =2{Ea + mA + ms)T{EA, rUs, itia) , 

92{q^) ={ms - mA)T{EA, m^, ttia) , (46) 
and where the function, T{EA,ms,mA), is, 



r^z? X - i{EA + ms)mA 



Thus, the ratio, 6'2(g^)/6'i(g^), gives. 



^i(g2) 2(EA + m, + mA) ' 
where, Ea-, the energy of A in the A^ rest-frame reads as 

1 



(4J 



EA = ^{ml^+ml-q') , (49) 
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with being the momentum transfer as previously defined. Working in terms of velocity, 
it may also be useful to express the invariant velocity transfer, tu(g^), as a function of the 
lambda masses, and nit^. Therefore, one obtains, 

Uj{q)=VA^-VA = , UJ.rain{q ) = I , (^max{q ) = ^ , (50) 

2mA(,mA 2mAi,mA 

where, va^ and v\, are the velocities of the baryon A;, and A, respectively. The effective 
QCD Lagrangian for heavy quark can be expanded both in powers of as{mb) and l/rrib- 
The radiative corrections will not be taken into account since they are not relevant in 
our analysis whereas the corrections proportional to AQco/fnb will be systematically 
calculated. These latter nonperturbative corrections are computed in the next section. In 
the following, all the form factors will be defined as a function of the invariant velocity 
transfer, u, instead of the momentum transfer, g^. 

4.2 The l/rrib corrections 

The QCD Lagrangian in HQET is written as 

Chqet = Kj^JvA, ■ Dhl^^ , (51) 

where, in the limit of infinite heavy quark mass, the quark field, h'^^ (x), replaces the 
heavy quark field, b{x), acting on x. Thus, the relation between the two quark fields 
reads as, 



hl^. i^) = exp 



irrih f At ■ X P+b{x) 



(52) 



with being the positive energy projection operator and Z)" being the gauge covariant 
derivative (= 9" — igst°'A'^). Including the corrections l/rnt to the effective Lagrangian 
previously written, one has: 



Chqet = K,. ivA, ■ Dh + 



2mft 



(53) 



where the gluon field strength is defined as usual by G"^^ = [iD°',iD'^] = igsf^G"^. In 
case of heavy to light quark mass transitions, the weak current will have the following 
general structure where the corrections l/rUb are taken into account: 

gT6^gT<^ + ^g.^<^+..., (54) 

where F reads as 7^ or 7^75. The dots define the higher order corrections which will be 
neglected in our analysis. 

By including the covariant derivative, D, as well as the corrections at the order of 
I /nib to the effective Lagrangian, it leads, respectively, to the local correction [T2] given 

by, 

SC^-' = ^mnl^ , (55) 
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and to the nonlocal corrections [Tz] given by, 

5C 
5C 
5C 



nlo,2 


- ' 






nlo,3 


- ' P 






nlo,A 


1 ~9 




2nib 2 



\2ub 



(56) 



where, g, holds for a light quark such as m, d and s, V is an arbitrary Dirac structure and, 
Qs is the field coupling constant. As usual, a^p is given by | [7a, 7/3]. 

4.2.1 The local correction 

Let us first start with the local term correction, to the effective Lagrangian, Chqet, 

written in the HQET formalism. It usually takes the form. 



(A(Pa, SA)\q TiDhl |Afe(pA,, saJ) = ua{pa, SA)(l)''{uj)T-f^UA,iPA„ sa 



(57) 



where one of the most general form of 

(f)''{uj) = (j)ui^^)vA, + 012(t^X + 013(^^)7 



U] IS, 



+ 



ifAt (p2l{uj)vl^ + (f)22{uj)v'^ + (t>23{uj)j^ , (58) 

because of spin symmetry properties. On one hand, the equation of motion for heavy 
quark. 



VA, ■ Dh 







applied to Eq. (jHTj) . gives 



VA, ■ {A{pA, SA)\q TiDhl^JAb{pA„ SA,)) =0 



(59) 
(60) 



and leads therefore to the following constraints where F = 1 and F = 75 have been used, 
respectively: 



«a(pa, Sa) [va, ■ 0(u;)]uA,(pAi,, saJ = , 
ma(pa, Sa) [va, ■ (j){uj)-f^]uASPAti saJ = . 
Thus, two relations between the 0jj(u;)'s can be obtained and read as, 

(j)u{uj) +UJ(f)i2{uj) = -023 (t^) , 
(j)2l{uj) +UJ(f)22{i^) = -013 (t^) • 

On the other hand, the momentum conservation also implies that, 

(A(pA, SA)\td^iq Thlj\Ab{pA„ SA,)) = 

{AipA, SA)\tDlqTh' |Ab(pA„ SA,)) + {AipA, SA)mDXjMPA„ SA,)) 



(61) 



(62) 



{mA, - mb)v^'^ - niAV^ \{AipA, SA)\qThl^JAb{pA„ sa,)) ■ (63) 
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By using the equation of motion for light quark, 

{lip - mg)q = , (64) 
then changing F into 7^r and contracting the index fi, one gets, 

(A(PA, SA)\q YTtDX^^ I Afe(pA„ saJ) = 

(mA, - mb)v^' - (mA - mg)v^yA{pA, SA)\q-f^'ThljAb{pA„ saJ) • (65) 

By choosing F = 1 and F = 75, Eq. gives, respectively, to 



u-1 

and, 
'uj + 1' 

OJ 



{uj) - 021 (C^)) + 



2w + 1' 



UJ 



> + 



ijj 



023 (t^) 



(mA, - m,){Fl{yj) + F°(^)) - [m^ - m,){Fl{yj) + ^F°(a;)) 



(0ii(u;) + 021 (t^)) 



2cj - 1 



013 (^^) 



4cu- 1 



023 (t^) 



(mA, - mO{G\{yj) - G\{yj)) + (mA - m,){G\{yj) + ^G'O(^)) . (67) 



In Eqs. (jUBj) and the form factors, -Ff (tu), and, G1{lJ), are the zeroth order form 
factors. The relations between them and the form factors, Qi{ij\ in the limit of m^ going 
to infinity, read as follows. 



and, 



Fl{uj) =2^2(c^) , 
Fl{yj) =0 , 



G\{u:) =2^2(0;) , 
Gl{u^) =0 . 



(68) 



(69) 



From Eqs. (jMj) and ()67|) . it is possible to express the 0jj(ci;)'s in terms of the zeroth order 
of the vector and axial form factors, F^{iJ) and Gliuj) and one has, 



01l(c^) = 

ijj{uj — 1) 
2(cj2 - 1) 



cj(cj + 1) 
2(cu2 - 1) 



-(m, - mAj(F°(a;) + F°(^)) - (mA - m,)(F°(^) + .^^"(a;)) 



+ 



(m,-mAj(G°(u;)-G?(c.^)) + (mA-m,)(G?(u;)+u;G°(cu)) 



(7^ - 1) 



0123(t^) , 



(70) 



15 



UJ - 1 

2(cj2 - 1) 



(C.-7) 



2(cu2 - 1) 

(mb-mAj(G?(cu)-G°(cu))-(mA-m,)(G?(cu)+a;G°H) 



- 1) 



+ 

0123(t^) 



(71) 



021 (c^ 



_ CJ(CJ+ 1) 

' ~ 2(cu2 - 1) 
cc;(a; — 1) 



K - ^Ajli'l'M + ^^2 M) + (^A - m,){F^{u) + ujF^{u)) 



+ 



2(cu2 - 1) 



+ 



^mlcc') , (72) 



and, 



(^22(^) 



+ 1 



2(cj2 - 1) 



(m, - mAj(F°(^) + F°(^)) + (mA - m,){F^{uj) + uF^{u)) 



2(cj2 - 1) 



(m, - mAj(G°(^) - G'?(^)) + (mA - m^){G\{u) + ^G°(^)) 

(1-7^) 



+ 



(^2 _ 1) 



0123(^) • (73) 



In Eqs. ()70II73|) . the assumption of (pi^ioj) ~ 023 (t^^) = 0123(1^) has been made for simphcity. 
It is now straightforward to derive the local corrections, 5FI°'^{ijj) and 5G'-°'^(co'), to the 
transition form factors, Fi{(jj) and Gi{uj) and one gets. 



5F["^\uj 



2mb 



6F^°'\u) = — 
nib 

1 

nib 



SF!,'^^\u) 



(j)u{uj) + (j)i2{uj)2{uj + 1) - 021 (tu) + 022(t^) 
2011(6^) + 2co'0i2(co') + 021(0^) + 022(ti^) 
011 (tu) + 02l(cU^ 



(74) 



and. 



SG[°'\lj 



2mh 



5G^'(cu) = — 
nib 



6G["'\uj) = — 



rrib 



4>u{uj) + (f)i2{uj){2u - 1) + 02l(cu) + 022(t^) 
20ii(cj) + 2uj(j)i2{uj) - 021 (c^) + 022(t^) 
012(u;) - 022(C^) 



Moreover, assuming that. 



6F[°'\u;) + 5F^°''(cu) + SFp\u) 



i/o,l 



2nib 



(75) 



(76) 
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it allows us to derive the expression of 0i23(<-o') and one obtains the following solution: 



>123l'^j = ^ 



-e + {rrib - mA,){F^{cu) + F°{u)) + (mA - mg){F^{u) + ujF^{uj)) 



, (77) 



where the form factor, 0i23(to'), is expressed in terms of the zeroth order form factors, 
Ff{(jj) and G°(a;), respectively. 

4.2.2 The non-local corrections 

Next is the first non local correction, 5C"'^"'^, given by. 



(78) 



which does not contribute because of the equation of motion, v\^^ ■ Dh^^ = 0, of the heavy 



quark, b. Thus, it directly leads to 5Ff"''^{uj) = 5Gf''''^{uj) = 0. The second non local 
correction, SC^^"'^, given by. 



(A(Pa,5a)K / rfl/T{g<^(0),<5/:'^'°'3(|/)}|A,(pA„SAj) 



2mb 



(A(pA,SA)|gr/^^^ (0)|A6(pa„saJ) , (79) 



and being a singlet under spin symmetry, leads to two form factors, Af (cj) and \f{uj), 
which only renormalize the zeroth order form factors, F^{uj) and G^{uj). Namely and 
independently of F, one has. 



G°(u;) = G°(u;) + 



2mh 



(80) 



For simplicity, it will be assumed that Af (tu) = Xf{uj) = Aj(cj). Moreover, these 
corrections can be safely neglected as they will only appear at the order of I /ml. 
Finally, by means of spin symmetry, the magnetic operator leads to the last non local 
correction, which reads as, 

(A(pA,«A)|i J rfyr|g-<^(0),|^^^<7„;,G"^/i^^(y)}|A,(pA„«Aj) = 



1 ' ~ ' ~-^l+^A.\.a/3. 



^^-UA{pA,SA)ipaf}{uj)T\^ ) S ^At (PA(, , SA J , (81) 



where. 



fapi^^)^ flli^^) + i>Abfl2{i^)\^aP + </'2l(t^) +?^A6</722(<^)]^a/3 , 



(82) 
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with the tensors, Sa/s and taf3, defined as, 

s-^ = l|^V-7V} and t"^ = ^{ 



(83) 



Assuming that (pii{uj) = ip2i{^^) and (puiuj) = ip22{^^), the corrections to the form factors 
Fi{q^) and Gi{q^) are given by. 



and, 



5F^'°'\uj) 



1 



Arrih 



1 



Arrih 



3v9ii(li;) + v3i2(co')(3 - 2uj) 
V9ii(u;) + V5i2(t^) 



4mfe 
1 

1 

4m/, 



3v3ii(u;) - (pi2{uj){?> + 2lj) 



15) 



Next and finally is the determination of the form factors, (fij{uj), coming from the magnetic 
operator corrections. By solving Eq. where the form factors, Fi{u) and Gi{uj), have 
been replaced by their expressions given in Eqs. ()87II88|1 . it yields to the expressions of 
the form factors, ipij{u), as a function of F^{u) and G^i^u), respectively. 

4.2.3 The form factors at the order of l/m^ 



After having derived in the HQET formalism the local, SC''°'^, and nonlocal, 
corrections to the effective Lagrangian, Chqet, the full form factors, Fi{u) and Gi{u!) 
read as, 



i=2 

4 



Gi{uj) = G'°(^) + SG\°'\uj) + 



nlo,j 



(86) 



i=2 



Explicitly, the expressions of the form factors Fi{uj) and Gi{u!) up to the corrections l/rrib 
can be written as follows, 



Fi{uj)=F^{uj) 



2mh 



011 (CU) + 20i2(t^)(l + 2uj) - 02l(t^) + 022(t^) 
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2mb 



V + 



■>21 



^22 



[uj) + u;(4</)i2(u;) + (fui^^)) 



F3H=F3V) + -^ 
and, 



012(t^) + 022(t^) - -(V?ll(t^) + <^12{UJ)) 



17) 



2mh 



(f)u{uj) + 2(f)i2{uj){2uJ - 1) + 021 (tu) + 022(t^) 



2mh 



40ii(t^) - 202i(t^) + 2(f)22{uj) + t^(40i2(u;) - i^uiuj)) 



+ - ^12M) 



nib 



(f)i2{uj) - 022(u;) - -{(pn{uj) - ifiui^^)) 



(88) 



The evolution of the form factors Fi{uj) and Gi{u!) as a function of the invariant velocity 
transfer, u, in case of A;, AV where V holds for the vectors p^, u and J/tp is plotted in 
Figs. 121 and EJ It can also be shown that the invariant velocity transfer, uj, to which the 
form factors Fi{uj) and Gi{u) have to be evaluated are equal to a; = 2.573 for — > Ap", 
UJ = 2.572 for A;, Aa; and lj = 1.856 for A^ ^J/ip, respectively. 



4.2.4 The integral overlap 

The calculation of the overlap between the initial and final baryon wave functions will be 
done by making use of the Drell-Yann approach where the current matrix elements are 
calculated in the impulse approximation. The two relations to leading order in momentum. 
Pa,, are. 



^2(0;) 



+ 



rriA + rriA^uJ rriA, + uiau 

G2{oo) ^ G^{uo) 
rriA + niA.uj 111^. + LuniA 



CsX{uj) , 
-Csl{uj) 



^9) 



where the overlap integral, X(co'), is, 

I{uj) = [ dxi(l>s{xiy<^>b{xi) =X(u;,mAj . 
Jo 



(90) 



The Clebsch-Gordan coefficient, Cg, suppresses the integral overlap since the s(ud) state 
is the only one which may contribute to the transition A;, to A. As a first approximation, 
the A may be seen as a superposition of various quark-diquark configurations. In our 
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model, one will follow the assumption of taking Cg equal to l/y/S. Then, let us expand 
the overlap function T{uj, m^J into inverse powers of the heavy baryon mass l/m\^, i.e. 





1 















ml 



where the leading term overlap integral, I^^\uj), is 



1 r 2^ 

exp |— K 



1^ Hi{R) 



and the next-to-leading order correction [12], X'^^)(ti;), reads as. 



(91) 



(92) 



J«(.) 
1 

+ 



UJ 



- - V2K(ab + ab) + Kab/(V2uj) 
b 



2V2b 



uoiuo + 1) HAk) 



1 

7! 



(a6 + ab) 



HAk) 



H[{R) 



uoiuj + 1) HAk) 



1 , uj-l HAk) 



H — a;6 



(93) 



with, Hi{x), a function defined as, 

/oo 
d^; (z + xY e~ 
■X 

where R and n are 



and H'i{x) 



2 ^cu3(cu + l)i^i(/t) 
difj(a;)/dx , (94) 



K = K 



UJ 



2u 



with K = \/2ab . 



(95) 



The parameters, a, a and 6,6 are defined in Section 4.3. Note also that the scaling 
functions, I^^\uj), and, T^^\uj), obey the normalization conditions = 1 and 

J(i)(l) = 0. From Eq. (jHHI), one can calculate the normalization condition, which gives 



Gi(l) = -c,J(l) + ( 



(96) 

at = 1. In a similar way of Eq. (j89p . another relation between the asymptotical form 
factors, 6^1(0;) and 6^2 (tu), is given in terms of the overlap integral, X(u;), between the A;, and 
A hadronic wave functions. Therefore, by working in the appropriate infinite momentum 
frame and evaluating the good current components the form factors, 6i{uj) and 62{^^), 
are obtained as the following, 

-Csl{uj), 



9i 

02 



2(^A + m, 
2(Ea + m. 



Csl{uj) . 



(97) 
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4.3 The baryon wave function 

Working within a quark-diquark approach, the baryon wave-function is proposed as 



^,(x,ki) = iVf$,(x)exp 
where the function, is defined as 

^i{x) = iVfv^V(l-a;)exp \-b^iM^{x - xoif 



{9t 



(99) 



with xoi being the peak of the baryon distribution function that represents the momentum 
fraction, x, carried by the heavy quark in the quark-diquark picture of the baryon. The 
parameter, bi, is related to the root of the average square of the transverse momentum of 
k^. The index, i, holds for a given baryon, namely A;, or A. We also impose the following 
condition of normalization for both of the functions, k^) and such as. 



/ / dxd2k^^,(x,kir^,(x,ki) = i 

Jo Jo 



dx ^(xY^(x) 



(100) 



The constants of normalization, A^* and N^, are directly given by making use of Eq. pOOj) . 
The parameters, hi and xa, read as. 



Xni = 1 , and ai = a -\ 



(101) 



where Mi is the baryon mass. From the average square of the transverse momentum of 
k_L given by. 



(ki) 



1 fOO 



JO 



d^k_L |^'i(x, k 



2 m2i,2 



one can effectively constrain the parameter, bi, since on gets, 

2&?(kl) = 1 . 



(102) 



(103) 



Regarding numerical values, we will follow Ref . ^3] '■ assuming that the root of the average 
transverse momentum, \/kT, is of the order of a few hundred MeV, one will choose an 



average of 500 MeV that gives bi = \pl GeV^^. The parameters, b and 6, are equal to 1.4 
GeV~^ and 0.1 GeV, respectively. The parameters, a and a, are equal to 1.0 GeV and 
0.3 GeV^, respectively. The baryon quark distribution, Di{x), takes the following form. 



Diix) 



dV ^.(x,ki)2 = - a:)Vexp[-26^M^(x - x,f 



(104) 



and is plotted in Fig. Elin case of A;, and A^. Finally, it will be assumed that the model 
of the baryon wave function of A;, is also valid for A^. 
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4.4 The weak decay amplitude 

In tree approximation, the effective interaction, Ti^-^-^, written as. 



^ 10 
^ i=l 



(105) 



gives the weak following amplitude, ^(Ai,A2)(^fe ~^ factorized into, 
^(,„,,)(A, ^ AV) = ^fvEy{A\-s^,il - 75)&|A.)(,,,,^) 



Ml(A,^AV)-MZ{A,^AV)\, (106) 



where the intermediate amplitudes, A4'^'^{Ah AV), are 



Mlf{A, AV) = V^^A 



T,P aT,P/ 



(107) 



The CKM matrix elements, Vj^(V^f„), read as VubV*g{VtbV^*) and VcbV*g{VtbV*g) , in case 
of Ab Ap, Ab Auo and A;, AJ/ip, respectively, fv and Gp are the decay constant 
of vector meson, V and Fermi constant. The amplitude, Ay^lai), expressed in term of 
Wilson coefficients, are the tree and penguin electroweak amplitudes which respect quark 
interactions in Af,-decay. They are listed in Section 4.6. Finally, the baryonic matrix 
element, {A\sjf,{l - 75)&|Af,) 



Bf^ ^ s), depending on the helicity state, (Aa, \\ 



reads as, 



B 



Py I + ruA 

y Ea + rUA 
_]_( Pv 



"(Aa.Av) 



Ciu;) + 2(2(00) 



(Aa.Av) 



^/2\EA + mA 

1 / Py 
y^yEA + niA 

C+(^) + 



C{u;) + C{u) 
C-ico)-CHco) 
Ciuj) 



p2 

^y 



EviEv + mA) 



(AA,Ay) = (|,0) , 
(AA,Ay) = (-|,-l) 

(Aa, Ay) = (i, 1) , 

(AA,Ay) = (-i,0) , 



(108) 



where, Py and Ey are the momentum and energy of V in the rest frame of Af,. The 
form factors C^(^^) = Ci(^) i C2(t^) are defined for convenience. The baryonic matrix 
element calculated in terms of the form factors, Ci{uj), can be converted in terms of the 
form factors, Fi{uj) and/or Gi{uj) as defined in Eq. 1)411) . As a function of Fi{uj), the form 
factors (i{u!) are written as follows. 



Ci(^) 

C2{i0) 



F2{UJ) 



Ma;) + F3(c.)fl + ^ 



(109) 
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and they are plotted in Fig. IHlas a function of a;, the invariant velocity transfer. Results 
are similar to those given in Ref. in case of — »• KJ /ip but they are different in case 
of Ab ^ Ap° or A5 ^ Ku. 

In our numerical analysis, the experimental values of the mass, my, and the width, 
Tvi of the vector meson, V , will be used to generate the vector mass resonance, m„. 
Therefore, the hadronic matrix element, {K{pi^,S!\)V{j)v,sv)\'Heff\^b{pk,,,SK^)) will be 
modified as follows: 



ml - my + iVymv 



X 



(my| J2|0)(A(pa, Sa)| Ji|A6(pa^, SaJ) 



(110) 



where Jj denote the quark currents and r„ refers to the radiative corrections in a^. Thus, 
the branching ratio, SS^\K\) V^A), will read as. 



myV 



V 



ml^ -ml + iVymy i^'^^^'v) 



KV)- 



XV) 



[III) 



Furthermore, when we calculate the branching ratio for A^ p°A, we should also take 
into account the — u mixing contribution since we are working to the first order of 
isospin violation. The application is straightforward and we obtain: 



^^{Ab ^ p°A) oc 



mpoTpO 



MUAb^Ap')- 



MiiAb-^Ap') 



+ 



MUAb^Au) 



n 



{sp - ml) + imj:^ 



(112) 



4.5 Operator Product Expansion 

The Operator Product Expansion (OPE) jT3 USl UZI is used to separate the calculation 
of a baryonic decay amplitude, into two distinct physical regimes. One is called hard or 
short- distance physics, represented by Wilson Coefficients and the other is called soft or 
long-distance physics. This part is described by Oi{p), and is derived by using a non- 
perturbative approach. The operators. On, entering from the Operator Product Expansion 
(OPE) to reproduce the weak interaction of quarks, can be understood as local operators 
which govern a given decay. They can be written, in a generic form, as. 

On = {qiTniqj){qk^n2qi) , (113) 
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where Tni denotes a combination of gamma matrices. They should respect the Dirac 
structure, the colour structure and the type of quark relevant for the decay being studied. 
Two kinds of topology contributing to the decay can be defined: there is the tree diagram 
of which the operators are Oi, O2 and the penguin diagram expressed by the operators O3 
to Oio- The operators related to these diagrams mentioned previously are the following, 



q' 



qi 



q' q' 

Oj = -g7/.(l - l5)b Y e<?'^Y(l + 75)g' , Os = -qal^.{l - 75)&/3 Y '^I'l'p^''^'^ + 75)?! 



O9 = T^Qli^i^ - l5)bY^q'^'^^i'^ - 75)g' , Oio 



2?"7m(1 - 75)fe/3 Y ^9'^^7''(1 - 75)gl • 

q' 

(114) 



In the above equations, a and /3 are the colour indices, denotes the quark electric 
charge and q', the quarks {u, d, c, s) which may contribute in the penguin loop. 

4.6 Wilson coefficients 

The Wilson coefficients [16 , Cj(/i), represent the physical contributions from scales higher 
than /i (of the order of 0{mh) in 6-quark decay) and since QCD has the property 
of asymptotic freedom, they can be calculated in perturbation theory. They include 
contributions of all heavy particles, and are calculated to the next-to-leading order (NLO) 
in such a way that one can get some corrections 0{as) from the leading- log-order (LO). 
By definition, C (/x) (we remove for convenience the index i) is given by |TS1 Ej , 



C{^i) = U{^i,Mw)C{Mw) , 
where f/(/i, Mw) describes the QCD evolution and reads as. 



w 



1 + ^ J 

4:71 



a, {My 



w 



w 



4:71 



J 



:ii5) 



:ii6) 



with J the matrix element including the leading order and the next-to-leading order 
corrections, d is the anomalous dimension. The final expression for C{fi) in the NLO, 



with U°{fi,Mw) = (a,(Mvi.)/a,(/i)) is 



C{^^) 



1 + 



47r 



-J 



w 



47T 



(117) 



where 5 is a constant term which depends on the factorization scheme. To be consistent, 
the matrix elements of the operators, Oi, should also be renormalized to the one-loop 
order. This results in the effective Wilson coefficients, which satisfy the constraint. 



aim,){Oiim,)) = C',{0,) 



tree 



:ii8i 
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where (O^)*'^^'^ are the matrix elements at the tree level. These matrix elements will be 
evaluated in the factorization approach. From Eq. (|118|) . the relations between C[ and Cj 
are [11111312011^1, 

C3 = C's — -Ps/3 , C4 = C4 + Ps , 

C5 = C*5 ~ -Ps/3 , Cg = Cg + , 

where 

P. = ^C2f^ + GK,/i,g2)^ , and P, = ^(3Ci + C2) f ^ + GK, /i, g^)^ ,(120) 



Stt V 9 V - ' 97r ' ' V 9 

with 



(^(mc, /i, g^) = 4 / (ix x(x — l)ln 



— x(l — 
7^ 



Here is the typical momentum transfer of the gluon or photon in the penguin diagrams 
and the expression of (^(mc/i, g^) can be found in Ref. [221 , Finally, the values for C[ 
are given in Table 1 where we have taken ^^(m^) = 0.112, aem{^b) = 1/132.2, nib = 4.9 
GeV, and rric = 1.35 GeV. 

4.6.1 Explicit Wilson coefRcient amplitudes for A5 AV 

Finally, in the following one lists the tree and penguin amplitudes which appear in the 

given transition: 

for the decay A;, ^ AJ/^, 

y4j/^(ai,a2) = ai , (121) 
^7/11,(03, ■ ■ ■ , Oio) = 03 + 05 + 07 + 09 ; 

(122) 

for the decay Ab Au, 
^2^^(01,02) =01 , (123) 

v^^^(a3, ■ ■ ■ , aw) = ^ (4(03 + 05) + 07 + 09 j ; 

(124) 

for the decay A;, Ap°, 
a/2^Jo(oi,02) = oi , (125) 

3 

a/2A^o(o3,--- ,Oio) = -(07 + 09) ; (126) 
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4.7 — LJ mixing scheme 

The direct CP violating asymmetry parameter, acp-, integrated over all the available 
range of energy of the vr+vr" invariant mass is found to be small for most of the non- 
leptonic exclusive 6-decays when either the factorization framework is applied. However, 
it appears that the asymmetry may be large in the vicinity of a given resonance, namely 
the uj meson in our case. To obtain a large signal for direct CP violation requires some 
mechanism to make both sin 5 and r large. We stress that — u mixing has the dual 
advantages that the strong phase difference is large (passing rapidly through 90° at the 
LO resonance) and well known ^Hl |211 123 12^1 12Z1 12H]- In the vector meson dominance 
model [^ni, the photon propagator is dressed by coupling to the vector mesons p° and u. 
In this regard, the — uj mixing mechanism |M] has been developed. Let A be the 
amplitude for the decay A;, — p^{uj)K — ^ vr+vr" A, then one has, 

A={k 7r-7r+|iJ^|Ab) + (A ti' ti+ \kb) , (127) 

with and being the Hamiltonians for the tree and penguin operators. We can 
define the relative magnitude and phases between these two contributions as follows, 

A= (A 7r-7r+|i7^|A;,)[l + re^V'^] , 

A={K 7r+7r-|iJ^|Ab)[l + re'^e-'^] , (128) 

where 5 and are strong and weak phases, respectively. The phase (f) arises from the 
appropriate combination of CKM matrix elements. In case of 6 — > ci or 6 — > s transitions, 
is given by = sj:g[{VtbV^^) / {VubV*^)] or sj:g[{VthV*^) / {VubV*^)], respectively. As a result, 
sin0 is equal to sin a (sin 7) for 6 c? (6 — > s) , with a (7) defined in the standard way [2]. 
Regarding the parameter, r, it represents the absolute value of the ratio of tree and 
penguin amplitudes: 

(pO(cu)A|/7^|A,) 



(pO(a;)A|/f^|A,) ' ^^^^^ 

With this mechanism, to first order in isospin violation, we have the following results 
when the invariant mass of tt+tt^ is near the uj resonance mass, 

(Avr-TT+li/^IA,) = ^n,.t. + ^t, , 

(A7r-7r+|/7^|A,) = -^Ilp^p^ + %p . (130) 

Here ty (V = p or uj) is the tree amplitude and pv the penguin amphtude for producing 
a vector meson, V, Qp is the coupling for tt+tt", lip^ is the effective p — uj mixing 

amplitude, and sy is from the inverse propagator of the vector meson V , sy = s — my + 
imvTv (with the invariant mass of the tt+tt" pair). We stress that the direct coupling 
UJ TT+TT^ is effectively absorbed into Hpt^ jH21 inHl IHll IHSl IHE], leading to the explicit s 
dependence of Hp^^. Making the expansion Iip^^{s) = Ylp^^m'^) + (s — m^)n|,^(m^), the 
p^ — UJ mixing parameters were determined in the fit of Gardner and O'Connell we 
will use 3fJe IlpUml) = -3500 ± 300 MeV^ ^pU^I) = -300 ± 300 MeV^ and 

H^^(m^) = 0.03 ± 0.04. In practice, the effect of the derivative term is negligible. 
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5 General inputs 



5.1 Source of CP violation: CKM matrix 



In most phenomenological applications, the widely used CKM matrix parametrization is 
the Wolfenstein parametrization jHHlEni- One of the main advantages in comparison with 
the standard parametrization [101 1^ is its easy analytical derivation at any order of A. a 
Four independent parameters, A, A, p and r], are usually used to describe the CKM matrix 
and each of these parameters can be (in) directly measured experimentally. By expanding 
each element of the CKM matrix as a power series in the parameter A = sin^^ = 0.2224 
{6c is the Gell-Mann-Levy-Cabibbo angle), and going beyond the leading order in terms 
of A in a perturbative expansion of the CKM matrix, it is found that the CKM matrix 
takes the following form (up to corrections of 0(A^)): 



Vr 



CKM 



where 



-A + lA^X' 
AXHl 



' 1 



p-if]) 



1 - 
-AX^ 



\X^ 



X 

\XHl+AA^] 



\AX\l-2{p + i^)) 



AX^{p-iri) \ 
AX^ 

i-iA^A^ ; 



p = p 



1 



2 



and r] = r][l 



2 



(131) 
(132) 



and T] plays the well-known role of the CP- violating phase in the Standard Model 
framework. From the CKM matrix, expressed in terms of the Wolfenstein parameters and 
constrained with several experimental data, we will take in our numerical applications, in 
case of 95% confidence level, 

0.076 < p< 0.380 and 0.280 < < 0.455 . (133) 

The values for A and A are assumed to be well determined experimentally: 

A = 0.2265 and A = 0.801 . (134) 



5.2 Input physical parameters 
5.2.1 Quark and hadron masses 

The constituent quark masses are used in order to calculate the electroweak form factor 
transitions between baryons and one has (in GeV), 

m^ = md = 0.350 , m, = 0.500 , mt = 4.900 , = 1.300 . (135) 

For hadron masses, we shall use the following values (in GeV): 

mAj, = 5.624 , mA = 1.115 , mj/^ = 3.096 , 

TTipO = 0.769 , = 0.782 . (136) 
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5.2.2 Form factors and decay constants 

The baryon heavy-to-light form factors, Fi{k'^) and Gi{k'^)^ depending on the inner 
structure of the hadrons have been calculated in Section 4. 

The decay constants for vector mesons, /y, do not suffer from uncertainties as large 
as those for form factors since they are well determined experimentally from leptonic and 
semi-leptonic decays. Let us first recall the usual definition for a vector meson, 

c(V^(<?)|gi7/.92|0) = /ymyey , (137) 

where my and ey are respectively the mass and polarization 4- vector of the vector meson, 
and c is a constant depending on the given meson: c = \/2 for the p° and uj and c = 1 
otherwise. Numerically, in our calculations, for the decay constants we take (in MeV), 

/,o= 209, /. = 187, /j/^ = 400. (138) 

Finally, for the total A;, decay width, rA(,(= 1/taJ, we use ta^ = 1.229 ± 0.080 ps. 



6 Physical results and simulations 

According to the preceding sections, all the essential parameters to perform precise 
simulations can be determined. We will particularly focuss on the A polarization, P'^, 
the PDM element, pJ^Q, of the vector-meson, the helicity and CP violating asymmetry 
parameters, and acp, of A^-decay and finally the branching ratios BTZ{Kb — > AV) 
where V holds for p^^ui and J/^j. 

6.1 Branching ratios 

By means of kinematical analysis developed in Sections 2 and 3 and the factorization 
procedure developed in Section 4, it allows us to compute the branching ratios of 
Af, 1 Af, h.uj and A^ — > KJ/il). The decay width of any process like A;, KV is 

given by the following formula 



r(A, ^ KV) = fEA±JI!±) -^J\J2 ^A„A. (A, - AV) "dQ , (139) 

where, in the A^ rest frame, the V^-momentum takes the following form, 

\/ ["^L - {'>T^v + mA)'^][ml - (my - mA)^] 

i^^>'l - ■ 

In Eq. ()139|1 . E/y and Q are, respectively, the energy of the A baryon in the A^ rest frame 
and the decay solid angle. The electroweak amplitude, Axj^^Xy{Ab — AV^), is given in 
Eq. ()l()(i|l . In order to take into account the non-factorizable term coming from the color 
octet contribution, calculations have been performed by keeping the effective number of 
color, N^-^f , to vary between the values 2 and 3.5. Branching ratios, BTZ, have been 
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calculated in case of A;, AV where V is J/tp,p^ and u and one obtains the following 
results: 

Bn{Ab AJ/ij) =(8.95 , 2.79 , 0.62 , 0.03) x 10"^ , 
Bn{Ah Ap°) =(1.62 , 1.89 , 2.16 , 2.39) x 10"^ , 
Bn{Ab Auj) =(22.3 , 4.75 , 0.19 , 0.64) x 10"^ , (141) 

for N^^f = 2,2.5,3.0 and 3.5, respectively. It is worth noticing that the experimental 
branching ratio 0, Bn^^^^^Ab AJ/ip) = (4.7 ± 2.1 ± 1.9) x 10"^ agrees with the 
theoretical predictions for 2.0 < N^^^ < 3.0. Regarding A;, Ap^ and A^ — > Auj, no 
conclusion can be drawn without any experimental data. 



6.2 Special case oi — lj mixing 

It is well known that the uj meson decays into tt+tt^ pair with a branching ratio of 2.2% 
and thus mixes with the p° meson by electromagnetic interaction. So, this mixing can 
be taken into account for the computation of the branching ratio BTZ{Ab Ap°) on one 
hand, and for comparison with the conjugate process A^ Ap° on the other hand. This 
comparison will allow us to check the direct CP violation in the sector of beauty baryons 
if a notable difference is found between A^-decay and its charge conjugated one. 

Computing the branching ratio of the charge conjugated process requires that only 
the CKM matrix elements, Vij, must be complex conjugated, while the intrinsic tree 
and penguin amplitudes are left unchanged. Owing to the important p° width, TpO = 
150 MeV/c^, by comparison to the uj one, F^^ = 8 MeV/c^, the branching ratios 
B7l{Ai,{Ai,) —>■ A(A)p°(ct;)), which matrix element is given by relation (112), is estimated 
by making use of Monte-Carlo techniques. The tt+tt" invariant squared mass, Sp, is 
generated according to a relativistic Breit-Wigner ^3] where the mass Mji and width 
are identified with Mp and Fp, respectively. For each generated invariant mass, Sp, the 
CP violating asymmetry parameter, acp{sp), between the two conjugated channels can 
be defined as, 

""^'^'^^ = BUiAb) + B-RiAb) ■ ^'''^ 

In Fig. 7, the CP violating asymmetry is plotted according to Sp. A remarkable effect is 
seen in the mass interval [720 — 820] Me\^/c^, which includes the u meson mass. Whatever 
is the value of N^^^ , one observes a strong enhancement of acp{sp) when the tt+tt^ 
invariant mass is in the vicinity of TSOMeV/c^. 

In Table El the numerous results obtained from the M-C simulations are gathered in 
case oi — uj mixing for A;, —>■ Ap^{u) — > Avr+vr" decay: mainly these are the mean 
branching ratio, BTZ, the asymmetry parameter, acp, computed in the whole range of 
TT+TT^ mass and the same asymmetry, acp{uj), limited to the u mass interval. Because 
oi p^ — uj mixing, the branching ratio BTZ{Ab Ap^{ijj)) weakly increases in comparison 
to that of BlZ{Ab — Ap°). We also notice that the acpiyj) for A;, — >• Att+tt" always 
reaches its maximum value when Sp is in the range of the omega mass. However, its 
maximum value strongly varies according to the effective number of color, N^^^ . Due to 
a lack of experimental data for the branching ratio BTZ{Ab Ap^{uj)), there is no way 
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of constraining N^^-' , unfortunately. In any case, the actual result can be seen as a clear 
signal of direct CP violation between beauty baryon and beauty anti-baryon. 

6.3 Helicity asymmetry and angular distribution 

6.3.1 Ah AV decay 

With respect to the definition of the helicity asymmetry parameter, a^^, given in Section 
3, its numerical value depends on the nature of the vector-meson V: 

a%{Ap''{u)) = 19.4% , a%{AJ/^) = 49.0% , 

which leads to a complete determination of the polar angular distribution of the A hyperon 
in the rest-frame. As far as azimuthal angle 0a is concerned and owing to the unknown 
value of the parameter p^''_ , the angle 0a is generated uniformly in the range [0, 27r]. In 
Fig.lHlare plotted the cos6'a and 0a distributions of the A-baryon in A;, rest-frame, in case 
of Afc —>■ Ap^{uj). Results are only shown in case of p° vector since they are very similar 
for J/ip. 

6.3.2 A — > p7z~ decay 

From preceding relations, both polar and azimuthal angular distributions of proton and vr 
can be obtained thanks to the complete determination of the A polarization, V^, and the 
PDM element p^2- Again, these values strongly depend on the nature of the vector-meson 
V coming from A;, decays. One obtains the following results: 

= -0.167 , pt, =0.25 (J/^) , 
= -0.21 , =0.31 (p°M) . 

It is worthy noticing that, in the framework of the factorization hypothesis previously 
exposed, the non-diagonal matrix element, pf2, is real, which makes easier the kinematics 
simulations. In Fig. El are also plotted the cos 6 and distributions of the proton in A 
rest-frame in case of A^ —>■ Ap^{Lo). Similar results have been obtained in case of J/ip. 

6.3.3 V i+i-, h+h- decays 

The angular distributions of lepton or pseudoscalar meson in the vector-meson rest-frame 
essentially depend on the normalized PDM element pj^. The latter is related to the 
probability of the vector-meson V to get an helicity value A2 = 0. Its numerical values 
are p^Q^ = 0.66 and pdg = 0.79, for J/ip and p°, respectively.. 

Despite the complicated form of W2{02,4>2) diplayed in Eq. the distribution 

of the azimuthal angle 02 will be uniform in the angular range [0, 27r] . Finally, in 
Fig. ^1 are shown the cos^^r and cos 6^ distributions for A;, — > Ap°(ct;) Att+tt^ and 
Af, — > AJ /tp — > Ap+p", respectively. 
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7 Time-odd observables 



Previously, it has been shown that, in an appropriate frame related to each intermediate 
resonance, the normal component of the resonance vector-polarization is a Time-odd 
observable and its measure could be a signal of TR violation. However, experimental 
determination of Pj^\i = A, V) is a difficult task and it requires high statistics data. 
Luckily, other possibilities to test TR symmetry exist. They have been developed by 
several authors and they rely on the search for Triple Product Correlations (TPC) among 
physical parameters represented by expressions like a, ■ {bj x Ck), where vectors di, bj and 
Cfe are either momentum or spin |¥3 EHl EZ] ■ It is worth recalling that the spin of any 
particle is not a direct measurable quantity and, only polarization with respect to a given 
direction can be determined. 

In the following, emphasis is put on TPC built from final kinematics variables measured 
in the A;, rest-frame. Needless to say that the detailed calculations performed in this work, 
both the (geometrical) angular distributions and the (dynamical) form factors deduced 
from HQET, represent important ingredients for precise computations of the TPC or 
their spectrum. In the recent years, L. Sehgal ^| and L. Wolfenstein ^\ analyzed 
very thoroughly experimental data coming from the reaction 7r+7r^e+e~, and they 

interpreted some kinematics asymmetry as T-odd effects. However, these effects do not 
necessarily indicate a TR violation, provided that complementary hypothesis (like non- 
conservation of CPT symmetry) have to be done. 

Inspired by the pioneering work of these authors, we give a special care to the analysis 
of some particular angles. First of all, we recall how the basis vectors of the A;, rest-frame 
(defined in Section 1) transform under TR: 

ez — > +ez , ex — > ~ex , ey = ez x ex — > -eV • (143) 

Azimuthal angular distributions like 0a in A^ rest-frame, (pp in A rest-frame and 4'i{4>h) in 
V rest-frame are directly computed according to the analytical expressions developed 
beforehand. For instance, let us consider the azimuthal angle 0a- Whatever is its 
distribution, the two parameters cos 0a and sin 0a are both even under TR and their 
proof is straightforward. According to the mathematical expressions: 

^ PA^ez , ^ ^ . , ^ ,^ 

u = — — - , COS0A = ey ■ M , sm0A = ■ (ey X u) , (144) 

\Pa X ez\ 

it is straightforward to see that cos 0a and sin 0a do not change sign under TR. Thus, in 
order to detect any TR asymmetry, we must look for other genuine triple products which 
can be associated to some specific angles. 

Let us define riA and ny as unit vectors which are normal respectively to A and 
vector-meson V decay planes. Their expressions are given by: 

PpXPn Pi+xpi- ^ Ph+^Ph- 
ua = 7^ ^ , nv = 7^ ^ , or ny = — ^— . (145) 

\Pp X p^\ \pi+ X pi- 1 \ph+ X ph- 1 

Those vectors are even under TR operation and their azimuthal angles, respectively 
0riA and (j)itv, will be referred to as 0(ni)- Performing the same calculations as for 0a, it 
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is shown that cos0(„.) and sin0(„.) are both odd under TR. Their analytical expressions 
are given by the following relations: 

Ui = ^ , cos0(„.) ^ey -Ui , sin0(„ ) = ez • (eV x Ui) , rij = riA , , (146) 

|ez X nil 

and, according to the transformation of each vector by TR, both cos0(„j) and sin0(„;) 
change sign under Time Reversal. An important question arises: what could be the 
order of magnitude of these TR violation processes? The only way to answer this 
question is the method of M-C simulations which provides us the interesting spectra 
of cos (f)(n,) and sin in order to cross-check the TR violation assumption. Exhaustive 
M-C studies were performed by using all the allowed values of the input parameters of our 
model and they reveal strong correlations of the T-odd observables with the azimuthal 
angular 0a distribution. 

(i) By adopting the simple hypothesis of a fiat distribution for da/d4)\, no dissymetries 
are observed in any of cos or sin distributions. Their mean vahics are very close 
to zero and the absolute value of the asymmetry factor, AS{X), defined by, 

AS{X) c ^dX - £ ^dX , (147) 

with X = cos0(„.) or sin0(„.) is less than 0.4%. This only corresponds to statistical 
fluctuations. 

(ii) Other solution which has not yet been exploited would be to assign a realistic 
distribution for the angle 0a. Indeed, the latter directly depends on the non-diagonal 
PDM element, p^''_ , which, a priori, is non-equal to zero. Departing from the relation 
given in Eq. (11), the 4>\ angular distribution can be inferred and its analytic expression 
is: 

da/d(l) oc 1 + ^aAs (3f^e(pi^_) cos - ^Jm(pi*_) sin 0^ , (148) 

where = 0a; a as = 49% for AJ/i/j and 19.4% for Ap°(a;). Normahzing 
d(T/d(f) provides a probability density function which must be positive in the whole range 
of variation: < < 27r. Furthermore, the (normalized) PDM elements could 
not exceed 1 in absolute value. Taking account of these kinematics constraints and 
adopting a conservative point of view, we suppose that ^e{p'^_) and '^m{p']^_) have 
similar contributions and the following choice is made: 3f?e(p^'L) = —'^m{p^_) = \/2/2. 
Interesting results concerning the asymmetry parameters, AS{X), are obtained in both 
the two channels with a sample of 10^ generated events for each one: 

Afo -> AJ/tp ; A5'(cos0,1a) = 4-3% , and A5'(sin0^^) = -5.5% , 
Ab Ap^{uj) ; AS{cos(l)fiJ = 2.4% , and A5(sin0^^) = -2.7% . 

Same calculations have been performed for cos 0^^ and sin 0^^ where ny = Up, fij/^. The 
corresponding asymmetries, AS{X), are close to zero: they are < 0.4% and compatible 
with statistical fluctuations. What could be the origin of these asymmetry discrepancies, 
despite the fact that the normal vectors n\ and nv respectively to A and V decay planes 
play similar role? 
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In our opinion, it is caused by the very difference of the azimuthal angular distributions 
(f)p in A rest-frame and 4'£{4>h) in V rest-frame. The former obeys to Parity violation in 
A-decay and it is dissymetric while the latter is flat because of Parity conservation so 
that it can be inferred that Parity violation in a given process like A hyperon decay 
could be the origin of Time-odd effect in a subsequent process. Other questions arise too: 
what are the true values of the real and imaginary parts of p+'L? Are the chosen values, 
9fte(p+'L) = — $^m(p+'_) = overestimated and be the real source of the observed 

T-odd effects? Whatever the physical origin of this novel effect is, all these interrogations 
are raised in the framework of the Standard Model because the main hypothesis underlying 
our calculations do not include any TR violation parameter. 

8 Conclusion 

We have studied the decay process A;, AV{1^) where the vector meson V is either 
J/ip^p^ or u;, with the inclusion oi — uj mixing when it was appropriate. When the 
invariant mass of the tt+tt^ pair is in the vicinity of the uj resonance, it is found that 
the CP violating parameter, a^p, reaches its maximum value. In our analysis, we have 
also investigated the branching ratios, BTZ, polarizations, P, and direct CP and helicity 
violating parameters, acp and uasi respectively, for the same channels. 

Thanks to an intensive use of the Jacob- Wick- Jackson helicity formalism, rigorous 
and detailed calculations of the A^-decays into one baryon and one vector meson have 
been carried out completely. This helicity formalism allows us to clearly separate 
the kinematical and dynamical contributions in the computation of the amplitude 
corresponding to A;, — ^ AV"(1~) decay. The cascade-type analysis^ is indeed very useful 
for analysing polarization properties and Time Reversal effects since the analysis of every 
decay in the decay chain can be performed in its respective rest frame. In order to apply 
our formalism, all the numerical simulations have been performed thanks to a Monte- 
Carlo method. We also dealt at length with the uncertainties coming from the input 
parameters. In particular, these include the Cabibbo-Kobayashi-Maskawa matrix element 
parameters, p and 77, the effective number of colors, N^^^ , coming from the factorisation 
scheme we have followed and the phenomenological model of the baryon wave functions 
for Af, and A. Moreover, the heavy quark effective theory has been applied in order to 
estimate the various form factors, F.i{q^) and Gilq"^), which usually describe dynamics of 
the electroweak transition between two baryons. Corrections at the order of 0{l/mi,) 
have been included when the form factors were computed. In the calculation of 6-baryon 
decays, we need the Wilson coefficients, C {mh) , for the tree and penguin operators at the 
scale mfc. We worked with the renormalization scheme independent Wilson coefficients. 
One of the major uncertainties is that the hadronic matrix elements for both tree and 
penguin operators involve nonperturbative QCD. We have worked in the factorisation 
approximation, with Nc treated as an effective parameter defined as N^^^ . Although one 
must have some doubts about factorisation, it has been pointed out that it may be quite 
reliable in energetic weak 6-decays. The recent experimental measurement of branching 

paper quoted in [501 and based on the same "cascade- type" approach has been released when our 
paper was in preparation. The authors using Monte-Carlo program and helicity analysis have analysed 
hyperon decays including lepton mass effects. 
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ratio for KJ/ip has confirmed the theoretical calculations we made for a value of 

Nl^f close to 2.5-3.5. 

As regards theoretical results for the branching ratios A;, — > AJ/t/j, A^ — > Ap° and 
Af, Au, we made a comparison with PDG for A5 AJ/i/j where an agreement is found 
(Bn^^^iAb AJ/iP) = 4.7 ± 2.1 ± 1.9 x 10"^ and BU'^'iAb AJ/^) = 3.1 x IQ-^). 
For Afo — > Ap° and A^ — > Aa;, the lack of experimental results does not allow us to 
draw conclusions. However, we made two theoretical branching ratio (in unit of 10~^) 
predictions, from 1.6 to 2.4 and from 0.64 to 2.23 for A;, — > Ap^ and A;, — > Auj, 
respectively. Wc have to keep in mind that, because of the difficulty in dealing with 
non factorizable effects associated with final state interactions, which are more complex 
for baryon decays than meson decays, these latter results are more dependent on the 
effective number of colors, N^^^ . In the special case of p° — a; mixing, the branching ratio 
Aft — > Ap^{uj) Att+tt", has a mean value around 2.15 x 10~^. The direct violating 
CP asymmetry parameter, acp, has also been calculated and a clear signal of direct CP 
violation between beauty baryon and beauty antibaryon has been observed. We have 
determined a range for the maximum of direct CP asymmetry, amaxi as a function of the 
parameter, N^^^ and the tt'^'tt" invariant mass. Moreover, the mean value of the direct 
CP violating asymmetry for Af, — > Ap°(a;) — > Att+tt", has been computed as well. 

By making use of the helicity asymmetry parameter, aAs-, determined for A5 — >■ AJ/ip, 
and Af, Ap^{uj), respectively, it allowed us to a complete determination of the polar 
angular distribution of the A hyperon in the A^ rest-frame. In fact, the knowledge of 
the A polarization, V^, which depends on the nature of the vector meson produced, in 
addition to that of the PDM elements, pfj", made available the polar and azimuthal angular 
distributions of the proton (and pion) in the A rest frame. In a similar way, the polar and 
azimuthal angular distributions of leptons and pseudo-scalar mesons in the vector meson 
rest-frame have been also computed. From weak decays analysis, one knows that vector- 
polarizations of outgoing resonances (or some of their components in appropriate frames) 
are T-odd observables. In our case, the normal components of and are T-odd, 
respectively and this might lead to a clear signal of TR violation. It is also well known 
that Triple Product Correlations (TPC) between momentum and spin can be used to put 
in evidence TR violation. However, not all TPC's can be exploited on the experimental 
side, because essentially of the difficulties to measure spin(s). In our work, we have shown 
that new and unsuspected observables can be measured: angular distributions of the 
normal vectors to the decay planes of the intermediate resonances in the Af,-frame. We 
found that the magnitude of their effects is directly related to the A;, polarization density 
matrix (PDM) and more precisely to the non-diagonal elements, p+'L, appearing in the 
interference terms of the decay amplitude. However, despite our ignorance of the real and 
imaginary parts of p^^_, sensible effects of Time Reversal asymmetry can be measured, 
provided the latter are not too small with respect to unity. 

In our opinion, new fields of research like direct CP violation and T-odd observables 
indicating a possible non-conservation of Time Reversal symmetry can be investigated in 
the sector of beauty baryons and especially the A^-bayons which can be highly produced 
in the future hadronic machine like LHC. In order to reach this aim, all uncertainties 
in our calculations still have to be decreased, for example non-factorizable effects have 
to be evaluated with more accuracy. Moreover, we strongly need more numerous and 
accurate experimental data in A^-decays, especially the A^ polarization. We expect that 
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our predictions should provide useful guidance for future investigations and urge our 
experimental colleagues to measure all the observables related to A;, baryon decays, if we 
want to understand direct CP violation and Time Reversal symmetry better. 

Acknowledgments 

The authors are indebted to their colleagues of the LHCB Clermont-Ferrand (France) team 
and to the Theory Group for their support and the interesting discussions regarding the 
very promising research subject of Time Reversal. One of us, (O.L.) would like to thank 
W. Roberts for stimulating discussions at the early stage of the paper. Another author, 
(Z.J. A.), is very indebted to J. -P. Blaizot, Director of the ECT* in Trento (Italy), for his 
kind hospitality at the Center and the very pleasant ambiance in which an important 
part of this work has been achieved. Finally, a recent graphical interface, so-called 
JaxoDraw jHIj has been used for drawing the diagrams in Fig. 2. 



35 



References 

[1] Z. J. A jaltouni, E. Conte, and O. Leitner, Phys. Lett. B614, 165 (2005), 
|hep-ph/l)412116| 

[2] Particle Data Group, S. Eidelman et al, Phys. Lett. B592, 1 (2004). 

[3] B. Andersson, The LUND Model, Cambridge Monographs on Particle Physics (1998). 

[4] M. Jacob and G. C. Wick, Annals of Physics 7, 404 (1959). 

[5] L. Ballentine, Quantum Mechanics: A Modern Development, Singapore World Sci. 
(1998). 

[6] J. Werle, Relativistic Theory of Reactions, North-Holland Publishing Company- 
AMSTERDAM , 419 (1966). 

[7] J. Jackson, Resonance Decays, High Energy Physics, Les Houches lectures, edited 
by C. DeWitt and M. Jacob (Gordon and Breach, New York, 1966 (1965). 

[8] J. G. Korner, M. Kramer, and D. Pirjol, Prog. Part. Nucl. Phys. 33, 787 (1994), 
,hep-ph/9406359. 

[9] F. Hussain and G. Thompson, (1994), |hep^ph /9502241| 

[10] M. Neubert, Phys. Kept. 245, 259 (1994), pli^h/9306320l 

[11] X.-H. Guo, T. Huang, and Z.-H. Li, Phys. Rev. D53, 4946 (1996), |hi^h/9706402 
[12] A. Datta, Phys. Lett. B349, 348 (1995), |hep^h/9411306[ 

[13] B. Konig, J. G. Korner, M. Kramer, and P. KroU, Phys. Rev. D56, 4282 (1997), 
|hep-ph/9701212 i 

[14] C.-S. Huang and H.-G. Yan, Phys. Rev. D59, 114022 (1999), |hep^ph/98113Q3 
[15] A. J. Buras, Lect. Notes Phys. 558, 65 (2000), |hi^h/9901409 
[16] A. J. Buras, (1998), |hep-ph/9806471( 

[17] G. Buchalla, A. J. Buras, and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996), 
fhep-ph/95 12380 

[18] N. G. Deshpande and X.-G. He, Phys. Rev. Lett. 74, 26 (1995), |hi^h/9408404 
[19] R. Fleischer, Z. Phys. C62, 81 (1994). 
[20] R. Fleischer, Z. Phys. C58, 483 (1993). 

[21] R. Fleischer, Int. J. Mod. Phys. A12, 2459 (1997), |hep^ph/9612446 

[22] G. Kramer, W . F. Palmer, and H. Simma, Nucl. Phys. B428, 77 (1994), 
|hep-ph/9402227| 

36 



[23] R. Enomoto and M. Tanabashi, Phys. Lett. B386, 413 (1996), hep-ph/96062T^ 

[24] S. Gardner, H. B. O'Connell, and A. W. Thomas, Phys. Rev. Lett. 80, 1834 (1998), 
rhep-ph/9705453, 

[25] X. H. Guo and A. W. Thomas, Phys. Rev. D58, 096013 (1998), |hi^^ /9805332l 



[26] X. H. Guo and A. W. Thomas, Phys. Rev. D61, 116009 (2000), |hep^ph/9907370^ 
[27] S. Gardner, (1998), |hi^h/9809479 



[28] S. Gardner, H. B. O'Connell, and A. W. Thomas, (1997), |hep^h/9707414[ 

[29] J. J. Sakurai, Currents and Mesons, University of Chicago Press (1969). 

[30] H. B. O'Connell, B. C. Pearce, A. W. Thomas, and A. G. Williams, Prog. Part. 
Nucl. Phys. 39, 201 (1997), hep-ph/9501251, 

[31] H. B. O'Connell, A. G. Wilhams, M. Bracco, and G. Krein, Phys. Lett. B370, 12 
(1996), hep-ph/95 10425, 

[32] H. B. O'Connell, B. C. Pearce, A. W. Thomas, and A. G. Williams, Phys. Lett. 



B336, 1 (1994), |hep^h/9405273| 

[33] K. Maltman, H. B. O'Connell, and A. G. Williams, Phys. Lett. B376, 19 (1996), 
|hep-ph/9601309| 

[34] H. B. O'Con nell, A. W. Thomas, and A. G. Wilhams, Nucl. Phys. A623, 559 (1997), 
|hep-ph/"970 3248 

[35] A. G. William s, H. B. O'Connell, and A. W. Thomas, Nucl. Phys. A629, 464c 
(1998),| hep-ph /9707253| 

[36] S. Gardner and H. B. O'Connell, Phys. Rev. D59, 076002 (1999), ,hep-ph/9809224l 

[37] S. Gardner and H. B. O'Connell, Phys. Rev. D57, 2716 (1998), |hi^h/97Q7385[ 

[38] L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). 

[39] L. Wolfenstein, Phys. Rev. Lett. 13, 562 (1964). 

[40] A. Hocker, H. Lacker, S. Laplace, and F. Le Diberder, Eur. Phys. J. C21, 225 (2001), 
[hep-ph/0104062, 

[41] L.-L. Chau and W.-Y. Keung, Phys. Rev. Lett. 53, 1802 (1984). 

[42] S. Pakvasa, S. P. Rosen, and S. F. Tuan, Phys. Rev. D42, 3746 (1990). 

[43] J. D. Jackson, Nuovo Cim. 34, 1644 (1964). 

[44] N. Byers, Spin and Density Matrix of decaying states , CERN 67 20, hep-ph/0412116[ 



[45] W. B ensalem, A. Datta, and D. London, Phys. Lett. B538, 309 (2002), 
|hep^h/0 205009| 



37 



[46] T. M. Aliev, V. Bashiry, and M. Savci, Eur. Phys. J. C38, 283 (2004), 
^ hep-ph/0409275 

[47] C. C. et al, Nucl.Phys.B(Proc Suppl) 115, 263 (2003), |hi^h/0210067 

[48] L. M. Sehgal and J. van Leusen, Phys. Rev. Lett. 83, 4933 (1999), |hi^ph/9 908426[ 

[49] L. Wolfenstein, Phys. Rev. Lett. 83, 911 (1999). 

[50] A. Kadeer, J. G. Korner, and U. Moosbrugger, (2005), |hep^ph/0511019| 

[51] D. Binosi and L. Theussl, Comput. Phys. Commun. 161, 76 (2004), |h^h/0309015[ 



38 



Figure captions 

• Fig. n Afe decay in its transversity frame and helicity frames for the A and vector 
meson V decays, respectively. 

• Fig. El Tree, color-commensurate, exchange and bow tie diagrams which are allowed 
in transition —>■ AMi. The squared dot represents the electroweak interaction in 
the Standard Model. 

• Fig. 121 Form factor distributions, Fi{uj), for A^ AV. 

• Fig. m Form factor distributions, Gi{u!), for Af, AV. 

• Fig. ElBaryon wave-function distributions for A^ and A. 

• Fig. ini Form factor distributions, Ci{uj), for Af, AV^. 

• Fig. 7 Branching ratio asymmetry, acp{uj), as a function of the vr+vr^ invariant mass 
in case of Af, — >■ Ap°(u;) — An'^ir' and for A'^^-^-^ = 3. 

• Fig. IHl cos 6'a and 0a distributions of the A-baryon in Af, rest-frame, in case of 
A,^Apyuj). 

• Fig. inicos^p and (pp distributions of the proton in A rest-frame in case of A^ — 
ApO(cu). 

• Fig. ^Icos^^ and cos^^ distributions in vector rest-frame, for Af, —>■ Ap^{uj) —>■ 
Avr+TT^ and Af, AJ/ip Afi~^fi~, respectively. 

Table captions 

• Tabled Vector-polarization under parity and TR operations. 

• Table 121 Wilson coefficients related to current-current tree and penguin operators. 

• Table El Branching ratio, BTZ, and direct CP violating asymmetry parameters, acp 
(mean value) and acp{uj) (maximum value), for Af, —>■ Ap^{uj). 
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Figure 1: Left-handed: decay in its transversity frame. Right-handed: hehcity frames 
for the A and vector meson V decays, respectively. 




Figure 2: From left-up to right-down, tree (T), color-commensurate (C), exchange (E) 
and bow tie (B) diagrams which are allowed in transition A;, AMi. The squared dot 
represents the electroweak interaction in the Standard Model. It can be shown that the 
contribution of each diagram to the total amplitude follows T ^ C ^ E ^ B. 
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Figure 3: Form factor distributions, Fi{u), for AV. The full, dashed and dot-dashed 
lines represent Fi(ci;), ^2(0;) and F^i^uj), respectively. 
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Figure 4: Form factor distributions, Gi{u), for A;, AV. The full, dashed and dot-dashed 
lines represent Gi{uj),G2{w) and 6^3(0;), respectively. 
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Figure 5: Baryon function distributions for and A. The full and dashed lines represent 
the Afo and A distributions, respectively. 
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Figure 6: Form factor distributions, C^iioS)-, for A^ — >• KV . The full and dashed lines 
represent C,\{yj) and C2('-^) respectively. 
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Figure 7: Branching ratio asymmetry, acp(<^), as a function of the tt+tt invariant mass 
in case of A;, Ap^{u) Att+tt" and for N^f^ = 3. 
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Figure 8: cos^a (left) and 0a (right) distributions of the A-baryon in A;, rest-frame, in 

dN 
d(t>A 



case of Ab Ap°(a;). The Y axis means (left) and j£ (right). 




400 — 



2 O O Q — 



1 750 — 



1 250 — 



750 — 



250 — 




Figure 9: cos6p (left) and (f)p (right) distributions of the proton in A rest-frame in case 
of Afe ^ Ap°{lJ). The Y axis means (left) and ^ (right). 
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Figure 10: cosOj^ (left) and cosO^ (right) distributions in vector rest-frame, for A;, 
A/)°(u;) Att+tt^ and A^ — > AJ/ip — >• A/x^/^", respectively. The Y axis means ^ 
(left) and (right). 
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el 


Odd 
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Even 


Even 


Pl 


Odd 


Even 


Pt 


Odd 


Even 


Pn 


E\'eii 


ODD 



Table 1: Vector-polarization under Parity and TR operations. 
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q^/ml = 0.3 




/ml = 0.5 




C[ 


-0.3125 




-0.3125 




C2 


+1.1502 
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-H 
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+3.839 X 10-4 
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-1.023 X 10-2 - 2.963 x 10" 


-H 


-1.017 X 10-2 - 9.938 X 10" 


-H 


^10 


+1.959 X 10-3 




+1.959 X 10-3 





Table 2: Wilson coefficients related to current-current tree and penguin operators. 



TV/// 


2.0 


2.5 


3.0 


3.5 


{BU) (10-^) 


2.10 


1.98 


2.15 


2.40 




6.4% 


1.04% 


-0.3% 


-0.2% 


acp{uj) 


96% 


15.0% 


-7.5% 


-3.3% 



Table 3: Branching ratio, BTZ, and direct CP violating asymmetry parameters, acp (mean 
value) and acp{oj) (maximum value), for — > Ap°(a;). 
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